Effects of a nonlinear bath at low temperatures 
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We use the numerical flow-equation renormalization method to study a nonlinear bath at low 
temperatures. The model of our nonlinear bath consists of a single two-level system coupled to a 
linear oscillator bath. The effects of this nonlinear bath are analyzed by coupling it to a spin, whose 
relaxational dynamics under the action of the bath is studied by calculating spin-spin correlation 
functions. As a first result, we derive flow equations for a general four-level system coupled to 
an oscillator bath, valid at low temperatures. We then treat the two-level system coupled to our 
nonlinear bath as a special case of the dissipative four-level system. We compare the effects of the 
nonlinear bath with those obtained using an effective linear bath, and study the differences between 
the two cases at low temperatures. 

PACS numbers: 03.65.Yz, 05.10.Cc, 03.67.-a 
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I. INTRODUCTION 

The linear bath of harmonic oscillators is the most 
common model system to describe a quantum environ- 
ment or bath, leading to dissipation and decoherenc e^ i- 1 
In this paper we examine an environment which is 
the simplest possible quantum-mechanical non-oscillator 
bath^: a single two-level system subject to an oscillator 
bath. The effects of this nonlinear bath are analyzed by 
coupling it to a spin, whose relaxational dynamics under 
the action of the bath is studied by calculating spin-spin 
correlation functions. 

Nonlinear baths have received a great deal of inter- 
est in the field of quantum computation, e.g., in the 
context of superconducting devices^ and for spin-based 
proposals2i2iS. 

There are two ways to define the border between sys- 
tem and bath in our total system. The first is to consider 
one two-level system to be the system and the other two- 
level system belonging to the bath. The bath consists 
then of a two-level system and the oscillators together 
and therefore it is a non-oscillator or nonlinear bath. 
Studying the problem this way is relevant, since, in many 
physical situations where the precise nature of the bath 
decohering a given system is unknown, it is modeled as 
a linear bath, with some given correlation function. It 
is, therefore, desirable to understand in more detail the 
kind and magnitude of possible errors introduced by such 
an approximation, in cases where the coupling cannot be 
assumed to be weak. 

The other way to define the system and bath is as 
follows: taken together, the two-level systems form a 
four-level system that is coupled to the linear oscilla- 
tor bath. Such systems of two two-level systems cou- 
pled to linear oscillator baths, representing a formidable 
problem in itself, have been studied in-depth using dif- 
ferent methods, such as functional-integral approachesiS, 
master equationsiii^, and numerical calculations within 
the quasi-adiabatic propagator path integral mcthod^'^. 
In this work we use a different approach: the numerical 
flow-equation renormalization method. 



The flow-equation method itself was introduced by 
Wegner^'' and by Glasek and Wilson^^. Using the flow- 
equation method one works in a Hamiltonian framework. 
The whole procedure is non-perturbative, as it relies on 
a unitary transformation and has an energy-scale sepa- 
ration built in. The extension of the flow-equation for- 
malism to dissipative systems was developed in Refs. IT^ 
ITtI and The effective Hamiltonian simplifles due 
to the fact that the coupling between bath and system 
disappears and equilibrium functions can be calculated 
easily. The method is non-perturbative, e.g., neither a 
Born nor Markoffian approximation is applied, which is 
what makes it attractive and new. Some drawbacks are 
that the fiow equations for the observables are not closed, 
and a special ansatz is needed. Here, we choose a linear 
ansatz, which works only for low temperatures smaller 
than the typical low-energy scale of the system Hamil- 
tonian. What is new in our work is that we apply the 
method to a general dissipative four-state system. The 
fiow equations are designed to yield correlation functions 
for non-ohmic baths as well, which, as we will see be- 
low, makes it possible to treat the linear bath approx- 
imation to our nonlinear problem on the same footing. 
In Refs. ilJi, ,2^ |^ and i25. subohmic, superohmic, 
and a peaked structured environment, were studied with 
fiow equations. Similarly, the linear-bath approximation 
leads to a new peaked structured bath spectral density, 
to our knowledge, never studied before. The numerical 
integration of the fiow equations allows us to study the 
nonlinear bath model by analyzing the spin-spin corre- 
lation functions. The numerical results are obtained at 
zero temperature. These results are then compared with 
those obtained from a linear bath with a non-ohmic bath 
spectral density. The bath spectral density is chosen in 
such a way that the linear bath leads to the same re- 
sults as the original nonlinear bath in the limit of small 
enough system-bath coupling strength. The bath and 
system correlation functions are both calculated using 
fiow equations. By carrying out the comparison, we are 
able to study the effects of the nonlinear bath at zero 
temperature. 
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The remainder of this work is organized as follows: 
First we describe the model in Sec. ^ Then, we write 
down the flow equations for the general dissipative four- 
level system. The flow equations for the nonlinear-bath 
system (Sec. IIII|) and its linear approximation (Sec. IIV|) 
are both special cases of these general flow equations. 
Finally, we discuss the numerical details in Sec. and 
compare the results of the two approaches in Sec. IVII 



II. THE NONLINEAR BATH MODEL 

As a model for a nonlinear bath we study a two-level 
system S coupled to a bath, which consists of another 
two-level system B coupling to a linear oscillator bath 
F: 



H = es<J 



S B 



E 



blbk 



(2.1) 



Here the parameters es and Ag serve to define any de- 
sired two- level system S. This system is coupled to B via 
af , with the coupling strength between S and B being 
given by g. The coupling constants are connected to 
the bath spectral density in the following way: 



J{uj) = n A^.(5(w — LOk) = aujfiuj/LUc 



(2.2) 



with the bath cutoff lUc, the systcm-bath coupling 
strength a and the cutoff function f{x). The oscilla- 
tions of af {t) at the frequency 2 A are noisy, due to the 
action of the harmonic oscillator bath F. The dissipative 
dynamics of S can be characterized in terms of several 
different quantities. Here we will analyze the decay of 
the equilibrium correlator (erf (t)crf (0)) by calculating 
its Fourier transform: 

Kiicj) ^ - / e-* (af (i)af (0)) . (2.3) 

Kf^{uj) is real- valued, for zero temperature it vanishes on 
the negative real axis and the integral over all frequencies 
gives 1. 



III. FLOW EQUATIONS FOR THE 
DISSIPATIVE FOUR-STATE SYSTEM 

We will write down the flow equations for the most 
general four-level system coupling to a linear oscillator 
bath. Therefore many other systems can be analyzed as 
special cases of this general framework. The Hamiltonian 
given in Eq. H2.1|l is a special case of: 

H^Y. ^"''^"/^ + E E S„^(6, + 6t ) 

a/3 k af3 

^afj{bk - bl) + y^CJfc 



blbk 



(3.1) 



The Tiap are the tensor products of the Pauli matrices 
= (Tq, ® cr/3 and the Greek indices are always summed 
from zero to three. (Tq is the unit matrix. By using 
the flow-equation technique we approximately diagonal- 
izc the model Hamiltonian given in Eq. (|3.1() by means 
of a unitary transformation: 



H{l)^U{l)HU{l)'^ . 



(3.2) 



Here, / is the flow parameter which is roughly equivalent 
to the square of the inverse energy scale which is being 
decoupled. The unitary transformation can be written in 
a differential form 



dH{l) 
dl 



where 



(3.3) 



(3.4) 



Using the anti-Hermitian generator ry in the canonical 
form T] — [Ho,H] leads to the fixed point H{1 oo) — 
Hq. We have investigated two different generators rj and 
?7, leading to two different fixed points, viz.. 



Ho 



k 



bW 



the oscillator bath, and 



a/3 



k 



UJk 



blbk : 



(3.5) 



(3.6) 



Q/3 



the total Hamiltonian without interaction terms. Note 
that we only have to derive flow equations with the gen- 
erator ry. The flow equations using rj are contained in the 
latter, due to the similar structure of the generators. 

The commutator [rj, H] contains coupling terms which 
are bilinear in the bosonic operators. We neglect these 
terms by truncating the Hamiltonian after linear bosonic 
terms. The bilinear terms are included by modifying the 
generator with an additional bilinear term as foUowsiLSi: 



V = [H„,H] +J2vk,: [hk + bl){b, -hi): . (3.7) 

kq 

The parameters rjkq are chosen in such a way that the 
bilinear terms are not generated. Because this cannot be 
done exactly we neglect terms which have the normal- 
ordered form: system operator times a bilinear bosonic 
operator. Here the normal order is deflned with respect 
to the non-interacting Hamiltonian. 

Comparing the general Hamiltonian given in Eq. I|3.1|) 
with Eq. H2.1|l specifies the initial conditions for the flow 
equations listed in Appendices 1X1 and IbI The flow equa- 
tions are nonlinear ordinary coupled differential equa- 
tions and they have the same structure as those obtained 
in Ref. 21 due to the similar truncation procedure. The 
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renormalization of the bath modes tOk vanishes in the 
thermodynamic hmit and can be neglectedi^. 

In order to calculate correlation functions the observ- 
ables have to be subjected to the same sequence of in- 
finitesimal transformations as the Hamiltonian. The ob- 
servable flow cannot be closed and therefore a linear 
ansatz has to be chosen which is only valid for temper- 
atures smaller than the typical low-energy scale of the 
four-level system coupling to the linear oscillator bath. 
The most general linear ansatz which is possible and used 
in the following is given by: 



o^Y. + E E /^f SaM&fc + bl) 



k 

-^Y.Y.-f^'^p^bk-bi). 

k af3 



The flow equations for the observables 
dO{l) 



dl 



Hi).o{i)] 



(3.8) 



(3.9) 



are closed according to the same normal ordering scheme 
as above and are listed in the Appendix [CI The correla- 
tion functions for I — > oo are defined by: 



m)o) 



tr{e-/5«o} 



(3.10) 



With the non-interacting fixed-point Hamiltonian _ffo 
such correlation functions can be calculated easily. 

Below we present the comparatively short result for 
the correlation functions for the fixed point Hq. They 
are found to be: 

{0{t)0) = + h°'h°' + h'^h'" + h'^h'i 



where 



+E[("fc + 

k 




Ck +nke 


^k 1 


(3.11 


ck = lA'iA' + 


^^f^ + 




ij ij 




+ -r/." + 






i j i j 

^ki^k 




- ^^T^T - 


- 




ij ij 

K^k 


(3.12 








i j i j 

^k^k ■ 





The Fourier transform of the correlation function is given 
by 

+ Y^\ink + 1)&{lo - ujk)ck + nkS{uj + ujk)ck ,(3.13) 

k 

where the Fourier transform is defined by 
1 /"^ 

(00)^^- dte-' {0{t)0) ■ (3.14) 



With these formulas different correlation functions can 
be calculated. The above formulas can also be used to 
obtain correlation functions for different operators, i.e., 
(O^^^(i)O^^^). The only difference in the result for Ck 

is that now a sum of products of the form Mi,^''°Vi^'"'* 
occurs and therefore the number of flow equations is in- 
creased. Below we put O = uf , i.e., analyze the decay 
of the equilibrium correlator (af (t)crf (0)). The correla- 
tion function is determined by the initial conditions of 
the flow equations. The limit of zero temperature T = 
can be obtained by putting nk = in the flow equations. 



IV. THE LINEAR-BATH APPROXIMATION 

In this section we study an approximation to the non- 
linear bath. We do this by replacing the composite bath 
B and F by a linear oscillator bath acting on the sys- 
tem S. What we obtain is a spin-boson system with a 
non-ohmic peaked bath spectral density 



(4.1) 



k k 

where the bath spectral density is given by 



(4.2) 



The linear-bath approximation should become exact in 
the weak-coupling limit g <C 1 according to Ref. |3 For 
higher coupling g we expect deviations to the nonlinear 
bath even in the thermodynamic limit — > oo due to, 
e.g., the non-separable structure^^ of our nonlinear bath 
model. 

We start our discussion by calculating the bath corre- 
lation function. The Fourier transform of the correlator 
of i? = gaf defines the "bath spectrum" (BB)^. For 
zero temperature, (BB)^ vanishes for negative frequen- 
cies. The Hamiltonian of the bath alone is given by 



H 



BF^^<y^+<J^Y.^'^(bk+bl)+J2^k:blbk: . (4.3) 

k k 



This is the initial Hamiltonian for the flow equations. For 
the purpose of calculating the bath correlation function, 
one has to solve the symmetric spin-boson problem. This 
cannot be done in a closed analytical waji^i^SiSi. An al- 
ternative method involves flow equations, which we ap ply 
again. We use the flow equation system used in Ref. |2l| 
with a generator of the form r]. The truncation scheme 
includes (as before) all coupling terms which are linear 
in the bosonic modes. The flow equation method is not 
restricted to any particular bath type, i.e., the bath spec- 
tral density of the linear bath can be chosen freely. 

The results for the bath correlation function (BB)^ 
are shown in Fig. ^ (BB)^ consists of one peak at the 
frequency w 2A, for small coupling a. For larger cou- 
pling it is shifted towards lower frequencies. Because the 
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FIG. 1: The (BB)^ correlation function for different values 
of a. The parameters are: 3=1 and A = 1. The height of 
the a = 0.005 peak is 5.2. 



bath F is ohmic, the position of the maximum of (BB)^ 
is given approximately by 2A{2A/ujc)'^"^^^~^"\ which 
is a result obtained within an adiabatic renormalization 
scheme^. In our results, the relation is best fulfilled for 
a/(27r) < 0.025. The low-frequency behavior predicted 
by the exact zero-temperature Shiba relation can only be 
found by solving the asymptotic flow equationsSi^ which 
is not done in this work (note the deviation from linearity 
for small w <C 1 in Fig. ^ e.g., graph for a/{2n) = 0.1. 
Thus, the bath spectral density shows a resonance at a 
characteristic frequency and behaves ohmic at small fre- 
quencies. This is why we expect the results to be com- 
parable to the ones obtained in Ref. But we stress 
that here, in contrast, the linear bath is modified with 
a second two-state system, compared to an additional 
oscillator. 

To obtain the final result for the linear-bath approx- 
imation we use the {BB)^ correlation function as bath 

spectrum J{lij) = tt {BB) ^ and solve the flow equations 
for the spin-boson system with this non-ohmic spectrum. 
The results are discussed in Section IVTI 

To close this section, we make a short note regard- 
ing higher-order correlation functions with respect to 
the fiow equations used here. Higher-order correlation 
functions are important to understand the behavior of a 
nonlinear bath. For small enough system-bath coupling 
strength, the nonlinear bath is identical to a linear one24. 
In this limit, the two-time bath correlation function plays 
a crucial role. For higher coupling strength, higher-order 
correlation functions become relevant. The linear bath 
has the special feature that higher-order correlation func- 
tions factorize into two-time correlation functions, since 
due to the Gaussian nature of the linear bath, the cumu- 
lant expansion breaks off at second order. The nonlinear 
case has generally no such feature, but we will see that 
the correlation functions obtained with the flow equa- 
tions factorize. The reason for this is the choice of a 



linear ansatz for the observable flow: 

aS) = h'^^iy, + y^f(l)cj^(bk + bl) . (4.4) 

In the limit Z ^ 00 we find for the symmetric spin-boson 
system 

aS^<^) = ^^f{l'^'^Mhk + bl), (4.5) 

and H{1 00) — Hq — J^k^f'^l^k- For / 00 the 
observables behave like linear oscillator variables and all 
higher-order correlation functions factorize. 



V. DETAILS OF THE NUMERICAL 
CALCULATION 

The flow equations are coupled ordinary differential 
equations. The integration method which turns out to 
be the most appropriate one is the fourth-order Runge- 
Kutta method with variable step size^S. The reason for 
this is that there are regions where a certain subset of 
the differential equations does not change much anymore 
and larger integration steps can be used, see discussion 
below. 

We have used a linear energy spacing Slo for the bath 
modes ujk and a sharp cutoff function f{x) = 0(1 — 
x). The results do not depend on these choices. The 
bath modes cuk are then given by Wfc = {k — l/2)Suj 
an d the coupling constants Afc are found to be Xk = 
(•\/a/7r-\/fc — l/2)5w, where Slo — ujc/N and N is the 
number of bath modes used for the numerics. For the 
four-level system the number of differential equations is 
GAN + 32 in the most general case, 5A^ + 7 for the non- 
linear bath with £5 = (see Appendix ID|| . lOA^ + 14 for 
the nonlinear bath with eg 7^ and 3A^ + 2 for the sym- 
metric spin boson system. We have used up to = 10"* 
bath oscillators. For the particular choice es = 0, the 
correlation function at discrete frequencies is given by 
Kf^{LOu) = {nl'^y/6uj, for the generator rj. All the h"'^ 
coefficients vanish for increasing flow I 00, as expected. 
For the other generator 77, we have to evaluate 

KUlo) - 5] |(0|(/j"3Eo3 + /ii3si3)|n)|2j(^ + Eo ~ E^) 

n 

+ 5](^")'l^0|I]ii|n)|2j(cj + So - S„ - oju) , (5.1) 

k,n 

where are the eigenvalues of Hs and |n) are the eigen- 
vectors of Hs for I ^ 00. In the parameter regime consid- 
ered below, it turns out that the equilibrium correlator is 
given by the same expression as for the generator rj with 
high accuracy. 

If we choose the generator ry leading to the diagonal 
fixed-point Hamiltonian ffo, which is scale-independent, 
then, since no asymptotic scale is present, the fiow equa- 
tions will first decouple the high-energy modes and then 
the low-energy modes. The same feature is found for 
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the correlation functions, see Fig. |5J They are first de- 
termined for high energies and the low-energy behavior 
is calculated last. The frequency regions of interest are 
centered on the peaks which appear at the shifted tran- 
sition frequencies of the four-level system. The spectral 
function around these resonances is determined by a sta- 
ble flow away from the asymptotic regime, which must 
be treated in a different wayi^t^^^. In practice, the flow 
equation can only be integrated up to I* w {6uj)~^ . We 
flnd that the integration can be stopped for much lower 
values; the system is decoupled with high accuracy for 
I > 20, see again Fig. |21 

For the other choice of the generator, 77, we find a differ- 
ent behavior of the flow. The flow equations first decou- 
ple the frequency regions around the peaks. The peaks 
contain gaps which are closed only slowly with increasing 
I, see Fig.O Therefore, in contrast to the generator rj the 
flow equations have to be integrated up to much higher 
values I « 10^. 

For both choices of the generator there are regions 
which are decoupled faster than others. Let us discuss 
the consequences for the numerics. The coupling con- 
stants of the Hamiltonian, X'^^ , and k'^^ decay while the 
flow parameter I is varied from zero to infinity. If we take 
a closer look at the flow equations, we observe that the 
differential equations for a certain mode k become trivial 
as soon as the X'^^ and are approximately zero, e.g. 
for the generator rj the higher modes are first decoupled. 
In the end, changes take place only in a low-frequency 
region, see Figs.|2|21 This is why, in practice, all the dif- 
ferential equations outside of this region can be replaced 
by trivial ones in the program. This accelerates the cal- 
culation significantly. 

The correlation function is given as a bilinear form in 
the /i^^ and J^^^. The coefficients ^"^^ and i/'^^ turn 
out to be zero at certain points. For the special case of 
the nonlinear bath with £5 = 0, only remains, which 
vanishes (intersects the frequency axis) at different points 
uj. This feature is contained in the flow equations and was 
observed in Ref. ^3 for the dissipative two-state system. 
The correlation function, bilinear in /ij,^, therefore has 
unphysical zeroes, which constitute a finite-size effect, 
only disappearing in the thermodynamic limit. 

There are different possibilities to address this prob- 
lem. In Ref. it was observed that for a certain value 
of the ffow parameter, the decoupled two-state system 
behaved like a dissipative harmonic oscillator. Further- 
more, it was shown that for the dissipative harmonic os- 
cillator there exists a conserved quantity which could be 
added at a certain value of / in order to complete the cor- 
relation function. Another possible way to deal with this 
problem is to integrate to a large value of I, such that the 
higher frequencies of the right peak are decoupled, as in 
Fig. 12 and to continue the integration in a restricted fre- 
quency range with a higher resolution of the bath modes. 
Finally, only a tiny gap is left which can then be closed 
by hand. 

Using the generator rj spikes appear around this singu- 
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FIG. 2: Behavior of the flow using the generator rj: Fourier 
transform of the equilibrium correlator for increasing flow pa- 
rameter I (solid line), compared with the final result of a 
Markoff approximation (dotted line). We see that the higher 
modes are decoupled flrst, i.e., the peak on the right hand is 
visible already for 1 = 7. At / = 20 the sum rule is fuIflUed 
with an accuracy of ~ 0.5%. This result is better than the 
Markoff approximation used to calculate the influence of the 
linear bath F on the system S + B, which should lead to good 
results for a/{2n) = 0.01. Naturally, the sum rule is not valid 
at the beginning of the flow. The parameters are: A = 1, 
As = 1.2, g = l, a/(27r) = 0.01, and N = 8000. 

larity, which are due to the finite number of bath modes 
N, i.e., we assume that the phenomena will disappear 
for increasing N. These spikes, generated via an ampli- 
fying effect, are not completely understood. Probably, 
they occur always if the change in the correlation func- 
tion happens too abruptly. This suggestion is motivated 
by the observation that the same feature appears, if the 
cutoff is chosen too close to the characteristic frequency 
of the system. There, we find that, due to the intrinsic 
properties of the differential-equation system, a spike ap- 
pears at the cutoff frequency. The system always needs 
a certain frequency range to spread out smoothly. 

The sum rules are fulfilled with an accuracy of ~ 0.5% 
using rj and ~ 0.1% using the generator 77, when not 
otherwise mentioned. 

The bath cutoff was chosen to he lUc — 5, when not 
otherwise mentioned. The results depend on the choice 
of the bath cutoff ujc- However, the qualitative behavior 
is independent of the cutoff and remains the same for 
different large-enough lOc- 

Flow equations give the best results for small coupling 
a. This does not mean that the approach is perturbative 
in the usual sense, but the neglected terms in the trun- 
cation procedure are smaller in the low-coupling limit. 
The flow equations for our specific system obtained with 
the generator rj were observed to give accurate results 
up to Q;/(27r) = 0.1 and the generator rj was used up to 
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FIG. 3: Behavior of the flow using the generator rj: Fourier 
transform of the equilibrium correlator for increasing flow pa- 
rameter / (solid line), compared with the final result of a 
Markoff approximation (dotted line). We see that the fre- 
quency regions around the peaks are decoupled at last. The 
sum rule is not valid at the beginning of the flow, but finally 
it is fulfilled with an accuracy better than ~ 0.1%. The pa- 
rameters are: A = 1, As — 1.2, g = 1, a/{2-K) = 0.01, and 
iV = 2000. 



a/(27r) = 0.08. 



VI. COMPARISON OF NONLINEAR WITH 
LINEAR BATH 

The relevant parameters in our model are eg, A5, A, 
the coupling strength g, and the coupling strength a to 
the linear bath F. We choose the time scale such that 
A = 1. Furthermore, the results discussed in the fol- 
lowing have been calculated for Ag = 1.2, close to 1, to 
keep the decay strong without being in resonance. We 
set the temperature to zero (T = 0) and vary the two 
parameters a and g. 

To begin our discussion, we note some generic features 
of the results obtained for the two approaches: linear and 
nonlinear bath. Since K^^ioj) is essentially the Fourier 
transform of the relaxational dynamics, it consists of 
several "Lorentzian like" peaks. Their number is con- 
strained to be less than the maximum number of 6 tran- 
sition frequencies for S + B m. the case of the nonlinear- 
bath system. In practice, degeneracies between transi- 
tion frequencies and selection rules reduce that number, 
e.g., to 2 for the unbiased case. More peaks, than the 
one transition frequency expected for a two level system, 
are observed in the case of the linear-bath approxima- 
tion. It turns out that the second peak is induced by the 
peaked bath spectral density. To summarize, we have 
found for the unbiased case two peaks for the linear and 
the nonlinear bath. One can associate two different deco- 
herence times with the two peaks. Both approaches show 
a similar behavior for increasing coupling where we ob- 



serve a larger separation between the two peaks. For the 
unbiased case es 7^ 0, we expect at least three peaks 
for weak coupling a. Additional possible zero- frequency 
delta-distribution contributions show up for the biased 
case, due to the definition of the correlation functions. 
Moreover, in contrast to the high-temperature results, 
no "pure" relaxation is seen in the biased case. 

Let us now focus on the unbiased case, i.e., es — 0. 
In the limit of weak coupling, g — > 0, all that remains 
is a broadened peak at the transition frequency 2A5 of 
system S alone. In that limit, the results for the two 
approaches coincide, as expected. This is visible in the 
topmost plots of Figs. O and El where the Fourier trans- 
form of the equilibrium correlator is shown for increas- 
ing system-bath coupHng g, while a/(27r) = 0.05,0.01 
is kept constant. We expect that for small g the non- 
linear and the linear-bath results fall together, i.e., in 
this limit a weak-coupling approximation"^^ or a Markoff 
approximation'^'^ are valid and one can therefore not dis- 
tinguish between linear and nonlinear bath. Further- 
more, the bath spectral density was chosen such that the 
two approaches agree for small enough coupling. 

With increasing g, the peaks are broadened and 
shifted, and additional peaks may appear (see Fig. 
Indeed, the most notable difference from a master equa- 
tion used for coupling of S and the nonlinear bath is the 
appearance of a second peak around the transition fre- 
quency 2A of the two-level fluctuator B. In this way, 
the power spectrum of the bath fluctuations appears in 
the short-time behavior of the correlator of the system S. 
This behavior cannot be captured by a master equation, 
where only one peak is present. Here, for the parame- 
ters chosen the linear and nonlinear bath agree nicely for 
g — 0.25, rather to our surprise, very well up to 5 = 0.5. 
For g > 0.5 the deviations become significant. Increasing 
g leads to a frequency shift and a change in the width of 
the "original" peak at 2A5. These changes are due to the 
change in eigenfrequencies and eigenvectors of the com- 
bined system S + B. If we keep g = 1 constant and com- 
pare the results for different increasing a, see Fig.01 the 
differences between linear and nonlinear bath depend on 
the ratio g/a. A small a/(27r) = 0.005 or a/(27r) = 0.01 
corresponds to a structured bath (Fig. and the re- 
sults show that the peak shape is similar for nonlinear 
and linear baths, but the peak position is different. We 
stress that for higher coupling, e.g., a/(27r) = 0.05 and 
a/(27r) = 0.1, the nonlinear bath acquires a new struc- 
ture and the shape of the peak deviates significantly from 
the linear bath. Thus, the qualitative differences between 
linear and nonlinear baths are smaller for a structured 
bath. 

Since the linear bath takes the full bath spectrum 
(BB)^ as input, this spectrum may show up in the result 
for the system correlator K^^{u)), as is indeed the case. 
Figures m and ini demonstrate that this effect is indepen- 
dent of the a values considered here, in contrast to the 
infinite temperature limilj^, where this effect was most 
pronounced for small 7, where the bath spectrum has a 



7 



I ^ 

^ ~ 

W 4 

S 2 

I 

° 1 

1 ^ 

'■B 0.5 



a- 

W 




1 

0.5 

0, 



a/(27i)=0.005 

: M 


1 1 , 
linear bath Jl 

nonlinear bath ,'i H 

^ l^A 


a/(27i)=0.01 




a/(27i)=0.025 


1 1 

/ \ 

^^^^^^/\,^^^^ 


a/(27i)=0.05 


/ ^ - 
1 \ 


a/(27i)=0.1 





Frequency CO 



FIG. 4: Fourier transform K^^{lS) of the equilibrium corre- 
lator of (5"f (i), for different structured nonlinear baths with 
a/(27r) = 0.005,0.01,0.025,0.05, and 0.1 from topmost to 
lowest graph. The values of the other parameters are: A = 
1, As = 1.2,3 = li S'lid es = 0. Nonlinear bath: solid line, 
linear bath: dashed line. Dotted line: Markoff approxima- 
tion. In Sec. [3 we discuss the singularities appearing in the 
correlators, e.g., in the lowest graph at frequency uj ~ 0.8 and 
the topmost graph. These two graphs were obtained using 
the generator rj, N = 5000, and I = 500. The rest of the 
graphs were calculated with the generator rj, N — 2000, and 
1 = 2- 10^ maximal. 



relatively sharp structure. To emphasize this point, wc 
note that even the asymmetric shape of the bath spectral 
density (BB)^ is mapped to the equihbrium correlator, 
see, e.g.. Fig. |S| compared with Fig. ^ This is an effect 
of the non-ohmic structure of the bath, which could not 
be observed for an ohmic linear bath. It can further only 
be observed when 5* and B are in resonance, i.e. for 
low frequencies the bath spectral density behaves ohmi- 
cally. The second peak, which appears for higher cou- 
pling g, is therefore determined by the bath correlation 
function {BB)^ with one peak in the vicinity of the sys- 
tem transition frequency and another given by the energy 
scale 2A5 = 2.4. For increasing system-bath coupling g, 
both peaks are shifted away from each other. Within 
the linear-bath approximation, the peak around 2A5 is 
shifted to higher frequencies to a lesser extent and the 
second peak around 2 A is shifted to lower frequencies, 
more pronounced compared to the nonlinear bath. 

The results of the linear and nonlinear bath agree well 
up to intermediate coupling strengths, as discussed pre- 
viously. Nevertheless, there are discrepancies: In par- 
ticular, the shifts are different. In contrast to infinite 
temperature, varying the parameters g and a leads to 
a visible shift of the peaks for the linear bath. Further- 
more, the peaks become wider and more asymmetric, but 
not as pronounced as the peak around 2A5 for infinite 
temperature. For higher values of a, the linear bath, 
in general, shows less structure than the actual nonlin- 
ear bath, however, the asymmetric shape of the peaks is 
mimicked to some extent. 



linear bath 
nonlinear bath 




Frequency O) 



FIG. 5: Fourier transform Kf^{uj) of the equilibrium corre- 
lator of af{t), for different values of the system-bath cou- 
pling g = 0.25, 0.5, 0.75, and 1 from topmost to lowest graph. 
The values of the other parameters are: A = 1, As = 
1.2,a/(27r) — 0.1, and es = 0. Solid line: nonlinear bath, 
= 2000, / = 8 ■ 10* maximal, and the generator r; was 
used. Dashed line: linear bath, = 2000, / = 1000, and the 
generator rj was used. 



Finally, to close this section let us turn to the biased 
case eg ^ 0. Figure [7| shows the equilibrium correlator 
of (jf (i) for two different values of the bias: 65 — 0.5 
and 65 = 1. For comparison, the Markoff approxima- 
tion is shown, which should lead to reasonable results for 
a/(27r) = 0.01, but the agreement between the two ap- 
proaches is worse than for the unbiased case. We note 
that the flow equations do not describe the third peak. 
What is the reason for this? Flow equations are equiv- 
alent to a unitary basis change. Therefore, we expect 
the same results for each unitarily equivalent Hamilto- 
nian. As soon as approximations are made, things do 
change, since, depending on the representation of the 
Hamiltonian, different terms have different signiflcance. 
Our approximation scheme leads to better results for the 
unbiased case. This was also seen in Ref. |^ for the spin- 
boson problem, where shifts of the bosonic modes were 
introduced and tuned to an optimal point in order to 
deal with this problem. We leave such an analysis for 
the dissipative four-state system for future work. 



VII. CONCLUSIONS 

The first main result of this work is the derivation of 
fiow equations for a general dissipative four-level system. 
With this framework, not only the correlation functions 
of our nonlinear bath model can be studied, but correla- 
tion functions for any four-state system coupled linearly 
to an oscillator bath. An additional key aspect is that 
there is no restriction for this bath to have an ohmic bath 
spectral density. 
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FIG. 



Fourier transform if^ 



of the equilibrium corre- 



lator of o-J'(t), for different values of the system-bath cou- 
pling g = 0.25, 0.5, 0.75, and 1 from topmost to lowest graph. 
The values of the other parameters are: A = 1, As = 
1.2, a/(27r) = 0.01, and es = 0. Solid line: nonlinear bath, 
N = 2000, 1 = 2- 10^ maximal, and the generator rj was used. 
Dashed line: linear bath, N — 2000, / = 1000, and the gener- 
ator rj was used. Dotted line (falls together almost perfectly 
with the solid line) : MarkofF approximation used to calculate 
the influence of the linear liatli F on the system .9 + B. 
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As an application of these general flow equations, a 
model of a nonlinear bath was discussed, consisting of a 
single two-level system subject to a linear oscillator bath. 
Its efi^ect on another two-level system at zero tempera- 
ture has been analyzed and compared with the effect of 
a linear oscillator bath, Many numerical results for vari- 
ous special cases have been obtained and discussed. For 
the flow equations of the nonlinear bath, two different 
generators have been used and the behavior of the corre- 
sponding flows has been analyzed. 

For small system-bath coupling the equilibrium corre- 
lator contains only one peak which can be obtained us- 
ing a Markoff approximation. In contrast the linear bath 
can describe the second peak appearing for larger cou- 
pling strengths. At least two frequencies are therefore 
present in the time evolution of the two-state system, 
and furthermore, the dccolica'encc is strongly increased if 
the system and bath are in resonance. As expected, the 
linear-bath approximation fails for the regime of large 
system-bath coupling. In that regime the linear bath 
may undoubtedly represent a good approximation to our 
actual nonlinear bath up to couplings on the order of 
half of the system energy scale. In the strong-coupling 
regime, on the order of the system energy scale, the agree- 
ment of the peak shapes is qualitatively better, when the 
bath spectrum has a strongly peaked structure. Here, 
discrepancies between the linear and the original nonlin- 
ear bath became clearly visible. Although we have only 
discussed the simplest example of a nonlinear bath, it is 
likely that the linear-bath approximation might lead to 
reasonable resiilts in the intermediate-coupling regime, 
also for other, possibly more sophisticated systems. 



2.5 

Frequency (0 



FIG. 7: Fourier transform K^^{uj) of the equilibrium corre- 
lator of (jf (i), for two different values of the bias es = 0.5 
topmost es = 1 lowest graph. The values of the other pa- 
rameters are: A = 1, As = 1.2, a/(27r) = 0.01, 5 = 1, 
ujc = 10, 1 = 100, and the generator rj was used. For es = 0.5, 
N = 2000 bath modes were used, the height of the delta peak 
at zero frequency is 0.229 and the sum rule leads to a value 
1.048. For es = 1, N = 1000 bath modes were used, the 
height of the delta peak at zero frequency is 0.544 and the 
sum rule leads to a value 1.149. Solid line: nonlinear bath, 
dotted line: Markoff approximation used to calculate the in- 
fluence of the linear bath F on the system S + B. 
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APPENDIX A: THE FLOW EQUATIONS FOR THE HAMILTONIAN 



Below we list the result for the dissipative four-state system discussed in Sec. IIIII A sum over i, j,l,p, s,r, and t 
from one to three is suppressed in the formulas. 



- - 4 2^ I - - (^Aj^^ Afe + \ + \ + ^k + f^k f^k + '^k'^k + '^k '^k + '*fe '^k J^k 

k 

+ (A^^'A^' + A°^AnM.'*"' + (A^"Al" + A^PAi^)6,H4° + (A'^^A^.y + A^^A^f" 6,,, + A'^'A^'e.y + A^^X'^e^u)^} 

In , 1 

- [(A^PA^' + A°PA°')A°'-e,ye,™ + (A^^'A^'e^y + A^°Aife,H + A^^y + A^'\'^e,u)>^^e,rn 

+ (AOP«;0' + A'Ti'nf)Kl-epije,rn + (A^P^'epy + A^"4^'e,,, + A'^^'^vU + nfe^)^^ e^^n] (2n, + 1) 

+ (A^^'A^' + A^PAHA^^^ep/n + (A^^Al" + A'P\'^)te^esu + (A^^A^pz,, + A^OAl"e,H + A^PA°'epz„ + A^"Ai°e,H)A^f } 

- = 4^ { [ - - >^TK" - - + A^o + .^^i - A«^«' - A^«i^)e.,„(2n, + 1)].^, 

k 

- [(A^"Al° + A'PX'P)Xl°esHe,rn + (A^^A^p,, + A^OA^^'e,;, + A'f A^p,, + A'^ \[°esu)Xl' e„n 

+ (A^"40 + A^PKln4°esi.e»™ + {A^P^epi, + A^^^e,;, + A'P^'e^,, + A^^'4°e,H)4^,™] (2nfe + 1) 

- (AOPAj^'ep,, + A^o.ci^"e,,„ + A^'p^'^pU + A^'^'^sinK' - (A^°4° + A^^A^Jf )A°"6,;„ - (A«f 4' + A^f 4)Ai:^"epy 
+ (A°PA5J'6py + A^°A^'6,,„ + A"PA"^py + A^nl°e,,„)4.^- + (A^^A^' + A'^Xl^K' e,,, + i^^'^k + A^^X[n4"esin} 



- 4 2^ - Aj. Aj. - A^. Afe - Kfe - Kj. Kj. + -(Aj, Aj. + Kj. jeip„ej/„ 

+ ^ [(4°Ai" + «rAl° - Af4" - Xr4")eun + + <^A« - A^V^' - A™4')e,;™] (2^^ + 

- [(A^f A^' + A°PAnArepye,™. + (A^^AI^'epy + A^"Ai^'6,,„ + A-^A^^y + A^-' A«'e.,„)A°'-e,,™ 
+ (A^^Ajf + A^^Ajf ) A™ esHEirn + (A°PA^fep/„j + A'*"A'j!^'esii + A*PA°'ep/,„ + A''™A^''es;i)A)J."eirn 

+ (A'^Pk^ + A ^ Kj. )Hk ^plj^jrm + (A'^^Kj, e^y + A ^H^egin + A ep/j + A £sln)l^t: ^jrm 

+ (A^o^io + A^P4P)4"e,He.™ + (A^P^epz™ + A'^n'^e^u + A'^^'^pim + A^™4°e,i,)4°e.™] (2nfe + 1) 

- (AOf 4' + A^P4)Af ep,„ + (A0p4epz, + A^"4e,i, + A^P4'ep,, + A^^"40e,i,)ArQt„e,™ 

- (A^°«;lo ^ 4P)A°™e,,„ - (A°P<'ep,™ + A^04'"e,,„ + A"P4'ep„„ + A«'"40e,i„)A°° 

+ (A^f A^.' + AOf AH'tf Mm - (A°PAf ep,, + A^"Ai^'e,,, + A'f A^p,, + A'^ X[hsH)4^eunejrm 

+ (A^OAi° + A^PAtP)4"e,,„ + (A°PA^'epz„ + A'^'XiY^esm + A"PA°'ep,™ + A--Ai°e,,„)4"} 
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dl 
dX9" 



dXV:° 



dl 



~dl 



(Al) 

q 

4(A«f AJ' + AO0A°'^epO-e,™ - 4(A0PA^p,, + A^OA^^e.j, + A'f A^p,, + A^^Ai°e,,OA"^e,™ 
= + 4(4° A'o + ^ A'^>ii„a;, + 2 ^ %,Af 

4(A«OAio + A^PX'^)A-\ueirn - 4(A°f Afe^y + A^OA'.^e,;, + A'f A^^y + A^^ AiOe,H)A'-^ei™ 
-wiAr + 4((4*A"' + KfA'^^im + (4° A'™ + KrA'°)eun)u>k + 2 J] m<,^r 

4{A'PXl' + AOPAO')A"'-ep,,-e,>„ - 4(A"f A^'e^,, + A^OA',^e,,„ + A"PA°'epy + A^^ A«'e,i„)AO'^e,v„ 
4(A^°A«' + A'PXi^)A^"'eMeirn - ^iA'^Xfepim + A'°Xi"'esH + A'f A°'ep,„ + A^'"Ai°e,iOA'^°ei™ 



"^fc _ ,.2 ..00 oV^„ ,.00 
dl 



Q 

-u,l4^ - 4(A2'A0' + Xi'Ai%ir,iVk - 2 5^,?,feK0" 

4(A"f 4' + A^f «^')A"'-epz,e,>„ - 4(A0p4ej,, + A^°4e,a + A'^K^Cpij + A'^ Ki°esH)A^'^ejrn 



^ = -u^lnf -AiXfA^' + X^A^^)eun^,-2Y,V,,Kf 



4(A«040 + A^^K^f )A'-"e,,,,£,,„ - 4(A"P«-4'ep;, + A'^n'^e^u + A'^^'^pi, + A^' 4'^su}A'' e^m 
- -wlnr - 4((A0^A"' + XrA'^^i^ + (A* A"" + XrA'yunju^k - 2^ry,,4'" 

4(A0f4' + A'P4)A^-epijejrm - 4{A^^Kl'epij + A'^k'^CsIu + A^^nfepn + A'^ n'^esii)A'''^e 
4(A''0k10 + A«f 4^')A™e,iiei™ - 4(A0f 4ep,„ + A«04™e«H + A^f 4'ep,„ + A'^Kfesu)A^''eirn 



'-jrm 



APPENDIX B: THE r)kq COEFFICIENTS 



%9 



2uJk 



+ Xl'Kfeun (So„) + XfK[°eun (S„o) + Xi'^^eun (So„) + >^^K^^eun (S„o) )a;fc 
+ (A 

+ A0'4'e«„ (Eo„) + XfKi^eun (S„o) + 4'eii„ (So„) + A^^K^^e«„ (S„o) )u 



9 

2a;„ 



+ (K°*A5j'eiim (2jnm) + l^-fj' X^k' ('iln {^nm) + '^q'^ ^ ^jlm {^nm) + l^'q" X!k (-iln (^nm) 



(Bl) 
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irn (^nrn) 

+ (AOPk« + A^PK|')<epyej™ (So„) + (A"f + A^^^e.^ + A'f + A^^^Oe^nXe,™ (So„) 

+ (A^"4° + A^P4f )4°e,,,e„„ (S„o) + (A"^ 4epo- + A^°4e,H + ^'^^^e^i, + A^^K^^e^zO^-'e.™ (S„o) 

+ A"PK^'ep,, + A^^ 

+ (A°f Ac"' + A^f K^O'^^epijej™ (So™) + (A°f ^Cp,, + A^^^e.^ + A'^K^'epi, + k'^ e^V^ e^m {^On) 

+ (A^04° + A'PK[n^l°esHeirn (S„o> + (A°^'<epz, + A^°4e,i. + A^^ ^e^;,- + A^^4°e,a)4^ei™ (Sno) " 

+ 4a;, [(A^f A|' + AO^ A^')A7epi,e,>„ (S„„) + (AO^ A^'e^,,- + A«°A^^e«;„ + A"PA^p,,- + A'^ Xi^e,in)X^^^€jrm {^nm) 

+ (A«°Af + A«PA5f)A™ {Ilnm) + (AOf Ajfepi™ + A'"Al!"e,H + A^^ A^' 

+ (A«f A|' + A'>PXf)\'l'epij€jrn (Son) + (A^^Af e^j,- + A^OA^^e.H + A'^Xfcpij + A'^ Xfe^JXp^m {^On) 

+ (A«OAi° + A^PAjf )A;°e,i,e„„ (S„o) + (A°^'Ai'epz,- + A^^A^^e^H + A^^ A^p^, + A^^Ai°e,H)A^^'ei™ (S„o) 

+ (A^^A^' + A"f Ai;')Arepi,ej™ (E„™) + (A°f A^'epz, + A^OA^-''e,i„ + A"^ A^^'ep,, + A'nfesin)Xl!^e,rm {^nm) 

+ (A^"a;,° + A'PX'P)Xresue„n (S„™) + (A°PA;'ep,™ + A^°A^'"e«K + A'^ A^'e^;™ + A^"A!,°e,H)A^fc"e„„ (E„„> 

+ (A^f A^' + A°f A°')A°'-epye,v„ (So„> + {A^^^y^e^ij + A^°A;/e«H + A*f A^'e^i,- + A'^ X'^ esH)^^ e^m (So„) 

+ (A«OaJ,o + A'PX'^)XfesHeirn (S„o> + (A^^'A^'ep^,- + A^^A^^e^n + A^^ A°'epi,- + A'^ X'^ esii)Xl^ etm (S„o) ] | • 



The expectation values of Sq,^ are defined by (Sa/3) = tr{S„^ exp(— /3if5B)}/tr{exp(— /JiJ^^)}- Here, -ffss is the 
Hamiltonian of the four-level system. 



APPENDIX C: THE FLOW EQUATIONS FOR THE OBSERVABLES 



- {A"P4' + A^^Kt)nl'e,i, - (A^°4° + A^f 4^')40e,H - (A0p4epy + A^'^e^u + A^^ ^e^y + A^^i^e.^O^ 
+ (A^^'A^' + A^PXt)vl^e,i, + (A^"Aio + A^^-A^^)^* + (A^^A^^y + A^OA^'^.h + A^^-A^p,, + A^^X'^c^uK] 

k 

+ \^^^ + 44' - AiVf ]e,,„(2n, + l))a;fe 

- [(A^^'A^' + A"f An//M.f.™ + (A"n^po- + A^°A^'6.H + A^n^po■ + A'^'X'^esu)^J^^e,rn 

+ (AOP4' + A^f 4>0'-ep,,e,™ + (A"f 4ep,, + A^o^^'e^n + A^^ 4'epz, + A'^ n'^e,u)i^l^ e^m] (2nfc + 1) 

- (A°f 4' + A^f 4)40ep,„ - (A^°4o + A^P4^')/ire.H - (A°P4epz„ + A^"4"e,,, + A^f 4'epi„ + A^"40e,H)4° 
+ (A^PA^' + AO^^A^Vfepi^ + (A"f A^'e^^^ + A^OAj^e.^, + A^f A^^^^ + A^'^X'^e^uK + (A^°A«' + A^f Ajf 
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/I f f \Oi,,ni \00,nO \m,,Oi -vnO.OO J>i,,ni „00,,nO „m,,Oi ^nO,,0 
— * \ 2 V ^'fc ~ M'k - ^k l^k ~ '^k f^k ~ l^k^k ~ ^k ^k ~ ^k ^k ~ l^k ^k 



- [(AOf Af ep,, + ^''^X'^esu + A^f AO'ep,, + A^^ AiOe^^OMfe'e.™ + (A^°AiOe,H + A«PA[''e,,iO/4°e»™ 

+ (A*040e,H + A^P^e.H);^^^.™ + (A^^^Mj" + ^'""^^sU + A^^^'^py + ^'^e^uVk' e„.n] (2«fe + 1) 

- (AOP«^'ep,, + A^°4^"e.,„ + A"P«^,y + A^^-4°e,,„)/x°^- - (A«"4° + A^^ K'^)^ife.,n - (A^^^^ + A^^'4)^^-'"6,,, 

+ (A°PA^'e^y + A«°Ajfe,,„ + A^^'A^^y + A^^Ai°e«;„);.°^- + (A^^'A^' + A°^'A°>re* + (A^^Ai" + A.^^X'^)vl\,i^] 

/) \ ^ \ ^ ( \00, ,nm \Om,,nO \nO,,Om \nm,, 00 ^00, .nm ^Om,,nO ^nO,.Om ^nm, .00 

- * 2^ \ 2\ ~ ^ ~ "^k i^k ~ '^k l^k ~ '^k l^k ~ '^fe '^/s ~ '*fc ^/s ~ '^fe '^fe ~ '*fc '^fe 

+ [(K^Vi" + 4"^!° - A^o^™ - Ar4°)e.;n + (4>fc' + 4Vfe - A5?Vfe"' - Ar^^fc°')e.i™] {2nu + l))a;fe 

- [(A^f A^' + AOf A°')Mrepij-e^rm + (A°PA^'epz, + A^°Ai^e,z„ + A"PA°'ej„, + A^^ A«'e,z„)Mfe'^ejrm 
+ (A^OAio + A'^\t),iresueirn + (A^PA^'ep,™ + A^o^^e.H + A^PA°'ep,„ + A^™AiOe,H)Mf e.™ 
+ (A°PkO' + A'^nf)vYepi,ejrm + (A°f K^'cpy + A^o^',^e,i„ + A"f ^o'epz, + A«-'K«'e,,„)i.O'^e,„„ 

+ (A«0 40 + A^P^Jf + (AO^ 4ep,™ + A^o^l^e.H + A*P4'ep,„ + k'^ esuVk e„n] (2nfe + 1) 

- (A°P4' + A^f ep,„ + (A"f 4ep,, + A^^^e,,, + A'P4'ep,, + e,u)^A' eu^e^rm 

- (A^"4o + A^f Atjf )Aire«z„ - (A0P4''ep,„ + A^OKl™e,z„ + A"f ^e^,™ + A-""Ki"e,i„)^°" 

+ (A«f A|' + A°f A°>f epz™ - (A°f A^pz, + A^OA^^e,,, + ^'"X^^e^ij + A«^Ai°e,K)4'^eit„e,y„ 

+ (A^"Ai" + A^f AiP)i/re.;n + (A^^A^'epf^ + A^^Af + A"^ A°'epi„ + A'^X'^esmH''} 



dl 



2^my° (CI) 

9 

^ = 2(4^/1°' + 4h^%lnC0, + 2 

- 4[(A«f A|' + AOf A^')/iO'-Mie,>n + (A^^'At'e^i,- + A^^A^^e^H + A'^A^p,,- + A'^ X[°esu)h'^ej 



dl 



dii""" 



dl 



,nO 



4[(A"f Af ep,,- + A^°A'/e«H + A^^ A^'e^^,- + A'^Xi^CsuW^eim + {A'°Xi° + A'PX'^)h^°esue, 
^ = 2((4'/i"' + K^'/iO0e«m + (4°/i''" + 4™/i'°h;n)'^fc + 2 5^%,Mr 

9 

- 4[(A^PA^' + A"f A°')/j"'-epz,e,™ + (A°f A^'cp,, + A^"Ai^e,,„ + A"^ A^^,, + A^ni°e,,„)^° 
+ (A^OAi" + A'PXinh"''eMeirn + (A°^'A^'epj„ + A^OAi"e,a + A^f A^pj„ + A-"Ai°e,a)/i''°eirn] 
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^^k° _ o\^„ „00 

-dT - -^Z.^^^^-^ 

Q 

q 

- 4[(A0f + A^o^lf e,;i + A'^nfe^n + k'^CsuW^ + (A^°40 + A'^n'^We^ueirn] 



- 4[iA^P4' + A'PKt)h"-epijejrm + {A"^" e^i, + A'^'^^esm + A'^^^^Cpij + A'^ K^;^esin)h°^e 
+ (A«°4° + A^^'K^f )/i™e«Hei™ + (A°f + A^^Kfe^n + A^^ ^'e^,™ + A'^Ki^esuWeim 



jrm 



APPENDIX D: A SPECIAL CASE, THE BIASED SYSTEM 



Here we present the flow equations for the specific case eg = 0, which was mainly used in this work. 



dl 

dA^ 
dl 

dA22 
dl 

rfA33 

dl 



k 

+ 8(-A0iAf 42 + A^^XfKf - A^^Xfnf + A^^Xl'^nf)] 

= J2Ai4'XfuJk + A^'Xl^Xl^ - A^^Af + A^'nl'^f - A^^nfnf) 

k 

= Y^AinfXfoj, + A'^^XfXf - A-^^XfXf + A^'nfnf - A^'nfnf) 

k 

= Y,^{-KfXfLOk - A'^XfXl^ + A^^Af A°3 - A^'Kfnf + A^^nf^f) 

k 

= Y^^i-^^Xfu:, - A^'XfXf + A^^Af Af - A^^^^^^ + A^^K^nf) 

k 



^ = -44^A°ic..+44^A2V. + 2^r;.,A°^ 

q 

- 4A0iA°iAf + 4A33^°1A^i + 4A^°A^^Xf - AA^^A^^Xf 

^ = -u:lXl'-AnfA''uj, + A4'A'^u;, + 2j2vk,Xf 

q 

+ 4A°ia33a°3 - 4A33A33Af - 4A^°A^°Xf + 'iA^^A'°Xf 
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- 4A0iA0iKf + 4A22A0iAt2i + 4A1°A334i - AA^^^^^^^f 

, ,2„21 ^\31a10, , _i_ /I \03 a22, , o „ ^21 



dl 



+ 4A01 A^^Kf - 4A22A2241 - 4AIOAI041 + 4A33AiOKf 



+ 4A"iAf /i^^ - 4A^^A^i/ii^ + 4Ai"A^^/i^" - 4A22Af /i^o 

= $:(-24V^^a;. - 4A0iAf 4^ + 4A33Ar4^) 

fc 

= ^(-2«-M^ia;. - 4A^°A^V^^ + 4A-A°Vr) • 



dl 
dl 



(Dl) 



The first equation for tiie coefficient A°° constitutes an energy renormalization. It is not necessary to take tfiis 
equation into account for tlie numerical integration of the correlation functions. The r/kq coefficients are given by 



Vkq 



20Jk 



{Xfnf + Xfnf (S33) - Xfnf (Soi) - Xf nf (Sio) 



+ {Xf4' (S22> + Xfnf (S33) - Xf4' (Soi) - Xl'nf (Sio) 



(D2) 



- 2iVq [{ - Kfxf (S33> - nfxf (E22) + Kfxf (Eoi) + nfxf (Eio) )uk 
+ ( - .°^Af (S33) ^ .fXf + <Af (Soi) + nfXf (Eio) )c 

I A 22 ^21^02 /V \ a22 ^21^21 /V \ a33 ^02^02 /V \ 1 a33 ^02 ^21 /v \\ 

+ 4c.,( - A°i«^2^f (Eoi) + A°Vf 41 (S22) + Aio«f «f (S33) - Ai"«2i4i ^^^o) 

+ A'^fKf (Eoi) - A^^^^l^l (E22) - A'^KfKf (E33) + A'^KfKf (Eio) ) 

+ 4a;,( - A^'XfXf (Eoi) + A^^Af Af (E33) + A^OAl^Af (E22) - A'^XfXf (E^o) 

- A'^XfXf (E22) + A^'XfXf (Eio) + A^^A^ A03 (Eqi) - A^^^^f A^^ (E33) ) 

+ 4u;g{ - A'^'XfXf (Eoi) + A'^'xfxf (E33) + A^^Af Af (E22) - A^'^XfXf (Eio) 

- A^^AO^Af + A^^Af Af (Eio) + A^^A^^A^^ (Eoi) - A^^Af A^^ (E33) ) | • 

For A > 0, and As > the only non-vanishing expectation values (E^/j) are 

(Eoo) = (Eii) = l 

(Eoi) = (Eio) 



A + As 



(E22) = — (E33) = 



^g^ + (A + As)2 
9 

^g^ + iA + Asf 
I 
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